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Abstract. In this paper we present a detailed study of the intra-cluster reddening material in the young open 
clusters located around I ~ 130° using colour-excess diagrams and two-colour diagrams. The study supports the 
universality of the extinction curves for A > A,/, whereas for shorter wavelengths the curve depends upon the 
value of the Rduster (total-to-selective absorption in the cluster region). The value of Rduster in the case of NGC 
654, NGC 869 and NGC 884 is found to be normal, whereas the value of Rduster in the cluster regions NGC 1502 
and IC 1805 indicates an anomalous reddening law in these regions. In the case of NGC 1502 the value of Rduster 
is found to be lower (2.57 ± 0.27) whereas in the case of IC 1805 it is higher (3.56 ± 0.29) than the normal value 
of 3.1. Although the intra-cluster material indicates a higher value of Rduster in the NGC 663 region, the error 
in the estimation of Rduster is too large to conclude anything. It is also found that the extinction process in the 
U band in the case of NGC 663 seems to be less efficient, whereas in the case of NGC 869 the process is more 
efficient. 
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1. Introduction 

Photometric studies are one of the most valuable and effi- 
cient tools for determining the physical properties of open 
clusters (e.g. distance, age, etc) and the interstellar matter 
(ISM) within the cluster as well as along the line of sight 
of the cluster. The accurate determination of the distance 
. of star clusters is crucial for a wide range of astronomical 
studies. The interstellar extinction and the ratio of the 
total-to-selective extinction R = Ay /E(B — V) towards 
the cluster are important quantities that must be accu- 
rately known to determine the distances photometrically. 

Observations of star clusters in the near-infrared 
bands, especially of young clusters which are still embed- 
ded in remnants of their parental clouds, are very use- 
ful to study the extinction behaviour in these clusters. 
Observations of interstellar extinction due to the intra- 
cluster matter yield fundamental information on the op- 
tical properties of the particles responsible for extinction. 
The analysis of extinction curves in 20 OB stellar associ- 
ations indicates that big complexes are obscured by the 
same type of interstellar matter (Kiszkurno et al. 1984). 
On the other hand Krelowski and Strobel (1983) com- 
pared the extinction law in two rich stellar aggregates, 
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namely Per OBI and Sco OB2, situated in different parts 
of the galactic disk. They found that the average extinc- 
tion curves of the two aggregates seem to differ substan- 
tially. The differences are large in the far-UV, which leads 
to the conclusion that an average extinction law cannot be 
applied to all associations. In an analysis of two different 
complexes of young stars (Sco-Ori region and Perseus re- 
gion) Krelowski and Strobel (1987) found different shapes 
of extinction curves, supporting the assumption that the 
obscuring material in these two young star complexes have 
different physical properties. Consequently, in the case of 
young open clusters, the average extinction curve should 
not be used to correct spectral or photometric data for 
interstellar extinction. 

There is much evidence in the literature that indicates 
significant variations in the properties of the interstellar 
extinction and these refer mainly to high values of the 
total-to-selective extinction ratio R as compared to the 
normal value, 3.1, for the galactic diffuse medium (see e.g. 
Pandey et al. 2000 and references therein). Whittet (1977) 
reported that the value of R in the galactic plane can 
be represented by a sinusoidal function of the form R = 
3.08 + 0.17 Sin(Z + 175°), which indicates a minimum value 
of R at I Kt 95°. However, Tapia et al. (1991) found a value 
of R = 2.42 ± 0.09 for the cluster NGC 1502 (I = 143.7°), 
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Table 1. The details of the clusters used in the study 



Cluster 


/ 


b 


Distance 


Log ag« 


E(B-V) 


Available 








(PC) 


(Yr) 




photometric data 


IC 1590 


123.13 


-06.24 


2940 


6.54 


0.32 


UBVIcJHK 


Be 62 


123.99 


01.10 


2513 


7.22 


0.85 


UBV 


NGC 436 


126.07 


-03.91 


2942 


7.78 


0.48 


UBVIc 


NGC 457 


126.56 


-04.35 


2796 


7.15 


0.48 


UBV I c 


NGC 581 


128.02 


-01.76 


2241 


7.13 


0.44 


UBVIc 


Tr 1 


128.22 


-01.14 


2520 


7.43 


0.63 


UBV 


NGC 637 


128.55 


01.70 


2372 


6.96 


0.70 


UBV 


NGC 654 


129.09 


-00.35 


2422 


7.08 


0.85 


UBVRJc 


NGC 663 


129.46 


-00.94 


2284 


7.13 


0.82 


UBVhJHK 


Be 7 


130.13 


00.37 


2570 


6.60 


0.80 


UBV 


NGC 869 


134.63 


-03.72 


2115 


7.10 


0.58 


UBVRJcJHK 


IC 1805 


134.74 


00.92 


2195 


6.67 


0.80 


UBVRJcJHK 


NGC 884 


135.08 


-03.60 


2487 


7.15 


0.58 


UBVRJcJHK 


NGC 1502 


143.65 


07.62 


900 


7.03 


0.74 


UBVhJHK 



which is significantly lower than the average galactic value 
of 3.1. 

The above discussions indicate that the interstellar ex- 
tinction shows a large range of variability from one line 
of sight to another. A precise knowledge of the spatial 
variability of interstellar extinction is important for the 
following reasons; 

(i) Since the extinction depends on the optical proper- 
ties of the dust grains, it can reveal information about 
the composition and size distribution of the grains. 
Consequently, variation of the extinction from one direc- 
tion to other may reveal the degree and nature of dust 
grain processing in the ISM (Fitzpatrick 1999). 

(ii) Since astronomical objects are viewed through in- 
terstellar dust, the wavelength dependence of extinction 
is required to remove the effects of dust obscuration from 
observed energy distributions. 

Uncertainties of extinction estimates limit the accu- 
racy of dereddened energy distributions. Such uncertain- 
ties might be acceptably small for very lightly reddened 
objects but can become important for modestly reddened 
objects. The intra-cluster extinction due to the remains of 
the star-forming molecular cloud decreases systematically 
with the age of the cluster (cf. Pandey et al. 1990). 

While analysing the UBVI C data of NGC 663 we found 
that the V/(U — B) colour-magnitude diagram (CMD) for 
NGC 663 cannot be explained by a normal value of the 
extinction ratio E(U - B)/E(B - V) = 0.72 (see Sect. 3). 
Some clusters in the direction of NGC 663 also appear to 
show abnormal extinction properties. This motivated us 
to study the interstellar extinction in detail towards the 
direction of the cluster NGC 663. 



2. Data 

The online available database at 

|http:/ /obs www.unige.ch /web da/ by Mermilliod 
(1995) provides an excellent compilation of the data on 
open clusters. We selected a range of 120° < I < 150° to 
study the intra-cluster extinction in open clusters. The 
online database lists 103 clusters in this longitude range. 
The parameters, such as distance, E(B — V), age etc, are 
available for only 54 clusters. Seventeen clusters in this 
direction are quite young (log age < 7.5) and are suitable 
to study the behaviour of intra-cluster ISM. UBV CCD 
photometry for nine young clusters of the above sample 
is given by Phelps and Janes (1994). A casual look at the 
(U — B)/(B — V) two-colour diagrams (TCDs) given by 
Phelps and Janes (1994, hereafter PJ94) indicates that 
the photometric data of 3 clusters, namely Be 62, NGC 
637 and NGC 663, do not satisfy the main-sequence (MS) 
curve, where apparent UV excess can be seen for stars 
having spectral type later than A. PJ1994 concluded 
that this may be due to UV excess of pre-main-sequence 
(PMS) stars and this part of the (U - B)/(B - V) TCD 
is not suitable for reddening estimation. We will discuss 
(U — B)/(B — V) TCDs in ensuing sections. The present 
analysis has been carried out using photometric data of 
14 clusters and details are given Table 1. 



3. Extinction 

One of the main characterstics of the diffuse interstellar 
matter in the Galaxy is its irregular structure. This makes 
it difficult to map the extinction as a function of galac- 
tic longitude and distance. Young open clusters also show 
differential reddening due to the remains of the associated 
parental clouds. The differential reddening decreases with 
the age of the clusters (cf. Pandey et al. 1990). However, 
it does not show any correlation with the location of the 
cluster in the galactic disk. 
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Thus the extinction in star clusters arises due to two 
distinct sources; (i) The general interstellar medium (ISM) 
in the foreground of the cluster, and (ii) the localized cloud 
associated with the cluster. While for the former compo- 
nent a value of R = 3.1 is well accepted (Wegner 1993, 
Lida et al. 1995, Winkler 1997), for the intra-cluster re- 
gions the value of R varies from 2.42 (Tapia et al. 1991) 
to 4.9 (Pandey et al. 2000 and references therein) or even 
higher depending upon the conditions occurring in the re- 
gion. 

3.1. Extinction Curve 

Extinction has often been analyzed using a two-colour nor- 
malizations of the form E(X — V)/E(B — V). In the present 
work the following methods are used to derive the ratio of 
E(\-V)/E(B -V). 

a) Colour excess diagrams (CEDs); method 'A' 

The colour excesses of the stars in the cluster region 
can be obtained by comparing the observed colours of 
the stars with their intrinsic colours derived from the 
MKK bi-dimensional spectral classification. For this pur- 
pose we used the data available in the online catalogue by 
Mermilliod (1995). When multiple data points are avail- 
able for a star, we used those selected by Mermilliod and 
given as MKS in the online catalogue. The colour excess 
in a colour index (A — V) is obtained from the relation 
E{\ - V) = (A - V) - (A - V) , where (A - V) is the 
intrinsic colour index and A represents the magnitude in 
the UBVRIJHK pass bands. Intrinsic colours are ob- 
tained from the MKK spectral type-luminosity class rela- 
tion given by Schmidt-Kaler (1982) for UBV, by Johnson 
(1966) for VRI converted to the Cousin's system using the 
relation given by Bessel (1979), and by Koornneef (1983) 
for VJHK. 

The colour excesses E(U - V),E(I - V),E{J - 
V),E(H — V) and E(K — V) are shown as a function 
of E(B - V) in Figs. 1, 2, 3 and 4. The stars having H a 
emission features and stars apparently lying away from 
the general distribution are not included in the analysis. 
A least-squares fit to the data is shown by a straight line 
which gives the ratio of E(\ — V)/E(B — V) for the stars 
in the cluster region. The slope of the line representing the 
ratio E(X — V)/E(B — V) along with the error is given in 
Table 2. In general, the least-square errors are quite large. 
The reason for the large errors is mainly the small sam- 
ple and the small range in the E(B — V). For comparison 
the colour excess ratios for the normal reddening law (cf. 
Johnson 1968, Dean et al. 1978) are also given in Table 2. 



b) Two colour diagrams (TCDs); method 'B' 

In most of the cases the MKK spectral classification is 
available only for a few stars of the cluster, which make 
the CEDs quite noisy (see e.g. the CEDs of NGC 654, 
NGC 663, NGC 869 and NGC 884). The TCDs of the 
form of (A — V) vs (B — V), where A is one of the wave- 
length of the broad band filters (R, I, J, H, K, L), provide 
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Fig. 1. E(U — V) vs E(B — V) colour-excess diagrams. 
Straight line shows a least-square fit to the data. 
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Fig. 2. E(I — V) vs E(B — V) colour-excess diagrams. 



an effective method for separating the influence of the nor- 
mal extinction produced by the diffuse interstellar medium 
from that of the abnormal extinction arising within re- 
gions having a peculiar distribution of dust sizes (cf. Chini 
and Wargau 1990, Pandey et al. 2000). On these diagrams 
the unreddened MS and the normal reddening vector are 
practically parallel. This makes these diagrams useless for 
determining the amount of reddening, but instead, very 
useful for detecting anomalies in the reddening law. Chini 
and Wargau (1990) and Pandey et al. (2000) used TCDs 
to study the anomalous extinction law in the clusters M16 
and NGC 3603 respectively. Figs. 5, 6, 7 and 8 show TCDs 
for the central region of the clusters. We used the data of 
the central region of the clusters to reduce the contami- 
nation due to field stars. The stars apparently lying away 
from the general distribution are not included in the anal- 
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Table 2. The colour excess ratios E(X - V)/E(B - V) obtained from the CEDs. 



Cluster 


E(U-V) 
E(B-V) 


E(I-V) 
E(B-V) 


E( J-V) 
E(B-V) 


E(H-V) 
E(B~V) 


E(K-V) 
E\B — V) 


NGC 654 






-2.20 ±0.26 


-2.47 ±0.26 


— 2.71 ± 0.27 


NGC 663 


2.12 ±0.32 




-2.28 ± 0.25 


-2.79 ± 0.30 


— 2.99 ± 0.30 


NGC 869 


1.90 ±0.13 




-1.20 ± 0.46 


-1.75 ±0.48 


-1.59 ±0.50 


NGC 1502 


1.74 ±0.18 




-1.87 ±0.16 


-2.04 ±0.23 


-2.17 ±0.22 


IC 1805 


1.76 ±0.08 


-1.34 ± 0.04 


-2.74 ±0.30 


-3.30 ± 0.37 


-3.66 ±0.38 


IC 1590 


1.56 ±0.30 


-1.33 ±0.17 
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Fig. 3. E(J-V),E(H-V), E{K-V) vs E(B-V) colour- 
excess diagrams for the clusters IC 1805, NGC 1502 and 
NGC 663. 

ysis. The slopes of the distribution, m c i uste r, are given in 
in Table 3. The slopes of the theoretical MS, m norma i, on 
the TCDs, obtained for the stellar models by Bertelli et 
al. (1994) are also given in Table 3. The errors associated 
with the slopes are significantly smaller than the errors 
obtained in the CEDs. The values of the (A - V)/(B - V) 
can be converted to the ratio E(X — V)/E(B — V) using 
the following approximate relation; 

E(X-V) _ m clustar l E(\-V) ] 
E(B-V) ~ m normal * l E (B-V)S normal - 

c) (U — B)/(B — V) two-colour diagram and colour- 
magnitude diagrams; method 'C 

In absence of spectroscopic observations, the (U — 
B)/(B-V) TCD and colour-magnitude diagrams (CMDs) 
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Fig. 4. Same as Fig. 3 but for the clusters NGC 869, NGC 
884 and NGC 654. 



are important tools to study the interstellar reddening to- 
wards the cluster as well as intra-cluster reddening. 

In the case of (U - B)/(B - V) TCD a MS curve (e.g. 
given by Schmidt-Kaler 1982) is shifted along a reddening 
vector given by the ratio of E(U - B) / E{B - V) = X 
until a match between the MS and stellar distribution is 
found. In this method the reddening vector X plays an im- 
portant role. Over the years the observational as well as 
theoretical studies have used ( or 'misused', according to 
Turner 1994) a universal form for the mean galactic red- 
dening law. However, theoretical as well as observational 
estimates for the reddening vector X show a range from 
0.62 to 0.80 (cf. Turner 1994). Recently DeGioia Eastwood 
et al. (2001) also preferred a value of X=0.64 in the case 
of Tr 14. The variability of X indicates the variations in 
the properties of the dust grains responsible for the extinc- 
tion. High values of X imply a dominance by dust grains of 
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Cluster 


(R-V) 
(B-V) 


(I-V) 
(B-V) 


(J-V) 
(B-V) 


(B-V) 
(B-V) 


(K-V) 
(B-V) 


NGC 654 


-0.61 ± 0.02 


-1.20 ± 0.04 


-2.16 ± 0.25 


-2.39 ± 0.23 


-2.50 ± 0.25 


NGC 663 


-0.50 ±0.02 


-1.13 ±0.03 


-2.33 ±0.30 


-2.73 ±0.35 


-3.07 ±0.43 


NGC 869 


-0.60 ±0.02 


-1.12 ± 0.02 


-1.90 ±0.10 


-2.41 ±0.10 


-2.53 ±0.11 


NGC 884 


-0.62 ± 0.02 


-1.06 ±0.02 


-1.96 ±0.08 


-2.44 ±0.08 


-2.59 ±0.10 


NGC 1502 






-1.92 ±0.24 


-2.04 ±0.24 


-2.08 ±0.24 


IC 1590 




-1.36 ±0.09 








IC 1805 




-0.98 ±0.04 


-2.12 ±0.15 


-2.58 ±0.17 


-2.73 ±0.18 


Normal 
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Fig. 5. (R — V) vs (B — V) two-colour diagram. The 
data point shown by open circles is not included in the 
least-squares fit (see text). 



small cross sections while the small values of X indicate a 
dominance of dust grains of larger cross section (cf. Turner 
1994). Turner (1994) raised a question about the so called 
reddening-free parameter Q [= (U — B) — 0.72(£> — V)]; 
how can Q be reddening-free when the reddening vector 
X is different from 0.72? 

As we have mentioned earlier, the distribution of stars 
in the V/(U - B) CMD of NGC 663 cannot be explained 
by a normal value of X (i.e. 0.72). The (U - B)/{B - V) 
TCD and the V/(B - V), V/(U - B) CMDs for the cen- 
tral region of the cluster NGC 663 are shown in Figs. 9 
and 10 respectively. The data have been taken using the 
105 cm Schmidt telescope of the Kiso Observatory, Japan 
(for details see Pandey et al. 2002). In Fig. 9, where the 
dashed curve shows the intrinsic MS by Schmidt-Kaler 
(1982) shifted along X — 0.72, we find a disagreement be- 
tween the observations and the MS at (B — V) ~ 0.90. 
PJ94 explained the disagreement due to the presence of 
pre- main- sequence (PMS) stars having UV excess. In Fig. 
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Fig. 6. (/ — V) vs (B — V) two-colour diagram. 



9 stars with V < 16 are shown by filled circles. The appar- 
ent distance modulus for NGC 663 is (m — My)=14.4 (cf. 
Pandey et al. 2002); the stars with V = 16 (i.e. My=1.6) 
will have mass ~ 2.5Mq. Since the cluster has an age 
of ~ 10 7 yr, the stars of ~ 2.5M© should have reached 
the MS and no longer be PMS stars. The TCD for these 
stars also does not support a normal value of X. We find 
that the MS shifted along a reddening vector of 0.60 and 
E(B — y)=0.68 explains the observations satisfactorily. 

The value of X in the NGC 663 cluster region can 
further be checked by comparing the theoretical zero-age- 
main-sequence (ZAMS) with the observed stellar distri- 
bution in the V/(B - V) and V/(U - B) CMDs. Once 
the reddening is known from the (U — B)/(B — V) CCD, 
the ZAMS is shifted to match the blue envelope of the ob- 
served stellar distribution in the V/(B — V) and V/(U — B) 
CMDs. The ZAMS fitting for E(B - V")=0.68 and appar- 
ent distance modulus (m— Mv)=14.4 is shown in Fig. 10a. 
Figure 10b shows V/(U - B) CMD where ZAMS, shifted 
for E(U — .B)=0.49 (corresponding to the normal redden- 
ing vector X = 0.72) and (m — My)=14.4, is shown by 
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Fig. 7. (J - V), (H - V), (K - V) vs (B - V) two-colour 
diagrams for the clusters IC 1805, NGC 1502 and NGC 
663. The data point shown by open circles are not included 
in the least-squares fit (see text). 
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Fig. 8. Same as Fig. 7 but for the clusters NGC 869, NGC 
884 and NGC 654. 
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Fig. 9. (U—B)/(B—V) two-colour diagram for the cluster 
NGC 663 . The continuous curve represents the MS shifted 
along E(U — B)/E(B — V) = 0.60, whereas dashed curve 
shows the MS shifted along a normal reddening vector. 
Star having V < 16.0 are shown by filled circles. 



a dashed line, which clearly shows disagreement with the 
distribution of the stars. The ZAMS for E{U - B)=0.40 
(corresponding to X = 0.60 and E(B — V) = 0.68) nicely 
fits the blue envelope of the distribution. This supports 
an anomalous value X=0.60 for the slope in the NGC 663 
cluster region. 

(U-B)/(B- V) TCDs and V/(B - V) CMDs were 
used to find out the value of X in all of the 14 clusters ex- 
amined in the present study. The (U - B)/(B - V) TCD 
for 3 clusters namely Be 62, NGC 436 and NGC 637 is 
shown in Fig. 11. The values of the reddening vector X 
obtained from the TCD/CMDs are given in Table 4. The 
uncertainty in the estimated value of X arises due to un- 
certainties in the intrinsic colours (i.e. ZAMS), uncertain- 
ties associated with the observations and also uncertain- 
ties associated with the visual fit of the ZAMS to the ob- 
servations. The typical total uncertainty in the reported 
values of X is estimated to be ~ 0.05. 



4. The Value of R 

Table 4 indicates that the clusters Be 62, NGC 637 and 
NGC 663 show a smaller value whereas the clusters NGC 
869 and NGC 436 show a higher value for X . The remain- 
ing 8 clusters show a normal value for X. The data given 
in Tables 2,3 and 4 are used to estimate the weighted 
(according to associated errors) mean value of the colour 
excess ratios E{\ — V)/E(B — V) and the ratios are given 
in Table 5. 



Pandey et al.: interstellar extinction in open clusters 



7 




B-V U-B 



Fig. 10. V/B - V (the left panel), V/U - B (the right 
panel) CMDs for the cluster NGC 663. Star having V < 
16.0 are shown by filled circles. The continuous curve in 
left panel represents the ZAMS fitting for E(B - y)=0.68 
and apparent distance modulus (m — Afy)=14.4. The con- 
tinuous and dashed curves in right panel shows ZAMS 
shifted for E(U - B)=0A0 and E(U - B)=0A9 respec- 
tively. For details see text. 

Table 4. The value of the reddening vector X = 
E{U - B)/E(B - V) obtained from the TCD/CMDs 



Cluster 


X 


IC 1590 


0.72 


Be 62 


0.60 


NGC 436 


0.84 


NGC 457 


0.72 


NGC 581 


0.72 


Tr 1 


0.72 


NGC 637 


0.53 


NGC 654 


0.72 


NGC 663 


0.60 


Be 7 


0.72 


NGC 869 


0.95 


IC 1805 


0.72 


NGC 884 


0.72 


NGC 1502 


0.76 



We define a parameter r which is the ratio of [E(X — 
V)/E(B — V)] c iuster (the ratio of colour excesses in the 
cluster region) and [E(X— V)/E{B— V)\ norma i (the ratio of 
colour excesses for the normal reddening law) . The values 
of E(X — V) /E(B — V) given in Table 5 are used to obtain 
the ratio r for the cluster region and resultant value of r 
is given in Table 6. For the clusters where the data from 



U band to K band are available,r is plotted as a function 
of A" 1 in Fig. 12. 

Fig. 12 indicates that the extinction in most of the 
cluster regions seems to be normal at A > A/, except 
for the clusters NGC 1502 (where the colour excess ra- 
tios for A > Xi are less than the normal one) and NGC 
1805 (where the colour excess ratio at A > Aj are higher 
than the normal one). Recent studies support a universal- 
ity of the extinction curves for A > Xj (see e.g. Cardelli 
et al. 1989, He et al. 1995). It is suggested that the nor- 
malization should be done using the E(V — K) instead 
of E{B - V) (Tapia et al. 1991) because the E(V - K) 
does not depend on properties like chemical composition, 
shape, structure, degree of alignment of interstellar dust 
(cf. Mathis 1990 and references therein). 

Cardelli et al. (1989) found that the mean R dependent 
extinction law can be represented by the following relation 

^=«A + ^ (1) 

where a\ and b\ can be obtained from the relations given 
by Cardelli et al. (1989). The above relation can be writ- 
ten, in terms of wgEijA , as 

ff*E$ = R(ax -l) + b x (2) 

The ratio of total-to-selective extinction towards the 
cluster direction 'Rduster' is derived using the eqn. (2). 
The value of ^gZ^p) — gi yen m Table 5 is used to 
estimate the value of Rduster and the results estimate for 
Rduster are given Table 6. The clusters that have a broad 
spectrum of data are discussed below. 

NGC 654 

Sagar and Yu (1989) concluded that at wavelengths 
greater than 5500 A, the extinction is normal. The pres- 
ence of unusually well aligned interstellar grains indicated 
by the polarization measurements seems to increase the 
extinction in the U and B bands slightly (Sagar and Yu 
1989). In the present work we find a rather normal extinc- 
tion law. The value of Rduster is ~ 2.97 ± 0.30(c) which, 
within the error, is close to the normal value of 3.1. 

NGC 663 

Using the colour excesses E(V — K) and E{B — V) 
Tapia et al. (1991) found weak evidence for an anomalous 
reddening law with a value of Rduster — 2.73±0.20, which 
is marginally lower than the normal value of 3.1. However, 
they felt that the scatter in their data is too large to con- 
clude about the value of Rduster- Yadav and Sagar (2001) 
reported values for E(X - V)/E{B - V) ((A > Xj)) which 
are significantly smaller than the normal ones. 

The (weighted) mean value of Rduster = 3.50±0.40(<r) 
suggests a marginally anomalous reddening law in the 
NGC 663 cluster region but in the opposite sense to that 
reported by Tapia et al. (1991) and Yadav and Sagar 
(2001). More near-IR data is needed to determine the 
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Table 5. Weighted mean value of the colour excess ratios E{\ — V)/E(B — V). 



E(U-V) E(R-V) E(I-V) E(J-V) E(H-V) E(K-V) 

E(B-V) E(B-V) E(B-V) E(B-V) E(B-V) E(B-V) 



NGC 654 


1.72 


± 


0.07 


-0.65 


± 


0.03 


-1.35 


± 


0.07 


-2.18 


± 


0.25 


-2.42 


± 


0.24 


-2.60 


± 


0.26 


NGC 663 


1.60 


± 


0.07 


-0.55 


± 


0.02 


-1.28 


± 


0.04 


-2.43 


± 


0.29 


-2.84 


± 


0.33 


-3.08 


± 


0.36 


NGC 869 


1.95 


± 


0.07 


-0.64 


± 


0.03 


-1.28 


± 


0.02 


-2.16 


± 


0.16 


-2.54 


± 


0.15 


-2.65 


± 


0.17 


NGC 884 


1.72 


± 


0.07 


-0.66 


± 


0.03 


-1.21 


± 


0.02 


-2.29 


± 


0.09 


-2.61 


± 


0.09 


-2.77 


± 


0.11 


NGC 1502 


1.76 


± 


0.07 














-1.97 


± 


0.19 


-2.14 


± 


0.24 


-2.20 


± 


0.24 


IC 1805 


1.72 


± 


0.07 








-1.25 


± 


0.05 


-2.55 


± 


0.22 


-2.86 


± 


0.23 


-3.07 


± 


0.24 


Normal 


1.72 
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Fig. 11. CCDs for the clusters Be 62, NGC 436 and NGC 637. The continuous curve represents the MS shifted along 
the reddening vector given in Table 4 and dashed curve shows the MS shifted along a normal reddening vector. 



Table 6. Mean value of the ratio r as a function of wavelength. 





2.90 






1.56 






1.25 






0.80 






0.61 






0.45 






Rclu 




r 


Cluster 


U 






Rc 






Ic 






J 






H 






K 












NGC 654 


1.00 


± 


0.04 


1.08 


± 


0.05 


1.08 


± 


0.06 


1.07 


± 


0.14 


0.98 


± 


0.09 


0.95 


± 


0.03 


2.97 


± 


0.30 


NGC 663 


0.93 


± 


0.04 


0.92 


± 


0.04 


1.02 


± 


0.03 


1.07 


± 


0.11 


1.09 


± 


0.11 


1.10 


± 


0.11 


3.50 


± 


0.40 


NGC 869 


1.19 


± 


0.02 


1.07 


± 


0.05 


1.02 


± 


0.02 


0.98 


± 


0.05 


1.02 


± 


0.05 


1.00 


± 


0.04 


3.04 


± 


0.20 


NGC 884 


1.00 


± 


0.04 


1.02 


± 


0.05 


0.97 


± 


0.02 


0.97 


± 


0.03 


0.99 


± 


0.02 


0.98 


± 


0.03 


3.19 


± 


0.12 


NGC 1502 


1.02 


± 


0.03 














0.85 


± 


0.05 


0.83 


± 


0.05 


0.81 


± 


0.04 


2.57 


± 


0.27 


IC 1805 


1.00 


± 


0.03 








1.00 


± 


0.02 


1.12 


± 


0.09 


1.14 


± 


0.08 


1.16 


± 


0.07 


3.56 


± 


0.29 



Rcluster in the NGC 663 cluster region. Here it is interest- 
ing to mention that the behaviour of the extinction curve 
towards UV also deviates from the normal one. Fig. 12 
indicates a lower value for the E(U — V)/E(B — V) ra- 
tio, whereas Yadav and Sagar (2001) reported a normal 
value for this ratio. They supplemented their data with 
the photometric spectral types which are based on the Q 
method, where they adopted E(U - B)/E(B - V)=0.72. 
Presumably a dominance of photometric spectral determi- 
nation forced the ratio of E(U — B)J E(B — V) to a normal 
value. 



NGC 869 and 884 (h and X Persei) 



From the extinction curve analysis Johnson (1965) 
found a value of R c i us ter = 3.0 in the NGC 869 and NGC 
884 cluster region. Tapia et al. (1984) also reported a nor- 
mal reddening law in the cluster region. Recently Yadav 
and Sagar (2001) found that the E(\-V)/E(B-V) ratios 
for A > A j are smaller than the normal ones. 

The colour excess diagrams (method 'A') of NGC 
869 indicates somewhat lower values for the ratios 
E(J - V)/E(B - V) and E(K - V)/E(B - V) but the 
errors are large, whereas in the case of NGC 884 the er- 
rors in the estimation of colour excess ratios are too large. 
The reason for the large errors is a small range in the 
E(B — V")(~ 0.3 mag). Because of the large errors we 
have not used the colour excess ratios obtained from the 
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Various studies have been carried out to estimate the 
value of R in the cluster region of IC 1805, but the re- 
sults are not conclusive. Johnson (1968), Ishida (1969) 
and Kwon and Lee (1983) reported an anomalous red- 
dening law in the cluster region with values of Rduster 
~ 5.7, ~ 3.8 and ~ 3.44 respectively. Kwon and Lee also 
reported a regional variation in the value of Rduster , with 
a maximum value of Rduster = 3.82 ± 0.15 for stars lo- 
cated in the outer region and the minimum of 3.06 ± 0.05 
for stars located in the central region. Sagar and Yu (1990) 
found that the interstellar extinction law in the direction 
of most of the cluster members is normal. The colour ex- 
cess ratios E(X - V)/E(B — V) for A > Xj obtained in 
the present work indicate an anomalous reddening law in 
the cluster region of IC 1805. The value of Rduster is es- 
timated as 3.56 ±0.29, which also indicates an anomalous 
reddening behaviour in the cluster region. 
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Fig. 12. The ratio r as a function of A -1 . 
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5. The A\/Ay curve 

Extinction has normally been analyzed using a two-colour 
normalization of the form E(X — V)/E(B — V). However, 
the true nature of the variability of observed extinction 
may be hidden by the choice of normalization. The quan- 
tity A\/Av reflects a more fundamental extinction be- 
haviour than the E(X - V)/E{B - V) (cf. Cardelli et al. 
1989). The average colour excess ratios given in Table 5 
can be used to estimate the quantity A\jAy in the fol- 
lowing manner, 



CEDs of NGC 869 in the further analysis. The TCDs 
(method 'B') indicate that the ratio of colour excesses 
E{X - V)/E{B - V) in both the clusters for A > Aj are 
perfectly normal. It is interesting to note that the CED for 
NGC 869 yields E{U - V)/E(B - V) = 1.90 ± 0.13. The 
V/(B-V) and V/(U — B) CMDs also seem to support the 
ratio of E(U-V)/E(B-V) = 1.9. In a recent study Keller 
et al. (2001) have adopted E(U - B)/E{B -V) — 0.72 to 
fit the ZAMS to the stellar distribution on the V/(U - B) 
CMD of the stars in the NGC 869 and NGC 884 clus- 
ter region. We find that only in the case of NGC 884 the 
V I (U — B) CMD supports a normal value for the redden- 
ing vector X in the cluster region. 

NGC 1502 

Tapia et al. (1991) found the ratio E(V - K)/E{B - 
V) = 2.20 corresponding to Rduster = 2.42 ± 0.09, which 
is significantly lower than the reddening value. Yadav 
and Sagar (2001) also found that the colour excess ratios 
E{X — V) /E(B — V) for A > A j are significantly smaller 
than the normal ones. The colour excess ratios for A > Xj 
obtained in the present study are in good agreement with 
those reported by Tapia et al. (1991) and Yadav and Sagar 
(2001). The value of Rduster is estimated as 2.57 ± 0.27 
which is in good agreement with that obtained by Tapia 
et al. (1991). 

IC 1805 



A\/Ay = [E(X - V)/E(B - V)Rduster] + 1 

where Rduster is taken from Table 6. In Fig 13a the nor- 
malized extinction in the form A\jAy is plotted against 
A" 1 for the clusters NGC 654, NGC 663, NGC 869 
and NGC 884 alongwith the average extinction law for 
Rduster = 3.1 given by Cardelli et al. (1989). Fig 13a in- 
dicates that the agreement between observations and the 
extinction law by Cardelli et al. (1989) is good barring the 
Au/Av values for NGC 869 and NGC 663. In the case of 
NGC 869 the ratio Ajj /Ay is higher than the normal one 
whereas in the case of NGC 663 it is lower than the normal 
one. The value of A v /Ay = 1.46 in the case of NGC 663 
supports Rduster ~ 3.5. Whereas in the case of NGC 869 
the colour excess ratios indicate a perfectly normal red- 
dening law. It seems somewhat strange that in the case 
of these two clusters the ISM behaves at Xu in a different 
way. 

The extinction law in the direction of two clusters IC 
1502 and IC 1805 is found to be anomalous. The nor- 
malized extinction A\/Ay for these two clusters along 
with the p Oph dark cloud (data taken from Martin and 
Whittet 1990) is plotted in Fig 13b. The effect of vary- 
ing Rduster on the shape of the extinction curves is quite 
apparent at the shorter wavelengths for different environ- 
ments of the star forming regions. 

As we have discussed the extinction in the direction 
of star clusters arises due to the general ISM in the fore- 
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ground of the cluster and also due to the cloud associ- 
ated with the cluster. Various studies using OB type sin- 
gle stars support a value of R ~ 3.1 for the general ISM 
(Wegner 1993, Lida et al. 1995, Winkler 1997). The min- 
imum reddening, E(B — V) min , towards the direction of 
the cluster is representative of reddening due to the fore- 
ground dust. The slopes of the distribution of stars having 
E(B - V) < E{B - V) m in on the TCDs can give infor- 
mation about the foreground reddening law. In the case 
of IC 1805 (and NGC 654), where E(B - V) mm is 0.65 
(0.77), we used stars having E(B - V) < 0.80 (0.85) to 
estimate the foreground reddening presuming that star 
having 0.65(0.77) < E(B - V) < 0.80(0.85) are not much 
affected by the anomalous reddening law in the cluster re- 
gion. The colour excess ratios E(J - V) / E(B - V), E(H - 
V)/E(B - V), and E(K - V)/E(B - V) obtained are 
-2.00 ± 0.32(-1.98 ± 0.28), -2.58 ± 0.35(-2.47 ± 0.27) 
and -2.80 ± 0.39(-2.82 ± 0.36) respectively, which sup- 
port a normal reddening law in front of the cluster IC 
1805. We further combined data of all the clusters hav- 
ing stars with reddening E(B — V) < 0.50. We feel that 
the limit of E{B — V) < 0.50 safely excludes the redden- 
ing due to intra-cluster matter as the smallest E(B — V) 
for the NGC 869 and NGC 884 is « 0.50 ( e.g. Uribe 
et al. 2002). A least-squares fit to 8 data points having 
0.01 < E(B — V) < 0.50 gives the colour excess ratio 
E(J-V)/E(B-V) = -2.23±0.32, E(H-V)/E(B-V) = 
-2.62 ± 0.22, and E(K - V)/E(B - V = -2.79 ± 0.26. 
These colour excess ratios also indicate a normal fore- 
ground reddening law towards the direction of the clusters 
used in the present study. 

Several studies have pointed to the apparent concen- 
tration of stars with high R- values in the vicinity of star 
forming regions. This effect has, for example, in the 77 
Carina nebula (Forte 1978, The and Groot 1983), M 16 
(Chini and Wargau 1990) and M 17 (Chini et al. 1980, 
Chini and Wargau 1998) and it may be presumed to be 
characterstic of many more HII regions ( Winkler 1997). 
Winkler (1997) compared the value of R obtained for 
hottest stars in the Galaxy (spectral type 08 or earlier), 
which can be considered as indicators of regions with re- 
cent star formation, and he confirms that in the majority 
of the cases the stars with large R indeed seems to be in 
the vicinity of star forming regions. 

On the basis of the above discussions we presume that 
the anomalous extinction law in the direction of cluster 
IC 1805 is due to the intra-cluster material. 



6. Conclusions 

In the present work we have carried out a detailed study of 
the intra-cluster material inside the young open clusters 
within a range of galactic longitude 123° < I < 144°. 
We used three methods (cf. Sect. 3.1) to derive E(X — 
V) /E{B—V) and the results agree with each other in a few 
percent. It is found that the behaviour of extinction varies 
from cluster to cluster. The main results are as follows; 
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values written along the curves represent the value of 
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1. The extinction curves at shorter wavelengths depend 
upon the Rduster while they converge for A > Aj. 

2. The extinction behaviour in the case of NGC 654, NGC 
869 and NGC 884 is found normal, whereas in the case 
of NGC 663 there is some tendency for a higher value 
of Rduster ■ It is interesting to mention that in the case 
of NGC 663 and NGC 869 the extinction at X v is found 
different from the normal one. In the case of NGC 663 
the extinction process in the U band seems to be less 
efficient, whereas in the case of NGC 869 the process 
is more efficient. 

3. The cluster regions of IC 1502 and IC 1805 show 
anomalous reddening laws. In the case of IC 1502 the 
Rduster ~ 2.57 ± 0.16 which is in agreement with the 
value (~ 2.49 ± 0.09) obtained by Tapia et al. (1991), 
whereas in the case of IC 1805 the present study yields 
Rduster ~ 3.53 ± 0.25 which is in agreement with the 
value (~ 3.44) given by Kwon and Lee (1983). 

More near-IR data is needed for the clusters NGC 663, 
NGC 1502 and IC 1805, where an anomalous value of 
Rduster is found, to study the reddening law in detail 
in the these clusters. 
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